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Abstract 

A coupled Volterra system is proposed. The model can be considered as one of the integrable 
discrete form of the coupled integrable KdV system which is a significant physical model. Many 
types of cnoidal waves, positons, negatons (solitons) and complexitons of the model are obtained 
by a simple rational expansion method of the Jacobi elliptic functions, trigonometric functions and 
hyperbolic functions. 



I. INTRODUCTION. 

The Volterra system [l| 

O'nt ~ 0,niO-n-l ^ ^Ln+l) = 0, (l) 

is one of the famous integrable differential-difference systems which has been apphed in 
various physical systems such as the network, statistical physics and biology [2i] . Some types 
of exact solutions of the model have been studied by many authors (say |3|). 

It is also interesting that the Volterra system is one of the simplest discrete form of the 
KdV equation. In fact, if we write 

an = 1 + 6\ ({n - 2t)S, -6H] + 0{5^) (2) 



3 

= l + 6Mx, T) + OiS^), (3) 
a„±i = l + 6\{x±6, r)+0{6^), (4) 

then becomes the well known KdV equation 

^{ur + Quu^ + u^^^)5^ + 0{S^) = 0. (5) 

Recently, some types of integrable coupled KdV system have been derived from some 
different physical fields including the atmospheric dynamics jj] , Bose- Einstein condensation 
Isl and two-wave modes in a shallow stratified liquid ol. A common special interesting case 

n n 

of j^-pj has the following form 

ut + Gavvx + 6uu^ + u^^^ = 0, (6) 
vt + Qvux + Quvn, + Vxxx = 0. (7) 

Some kinds of analytical negatons, positons and complexitons of the coupled KdV system 
©-([Tj) for a = —1 are studied in [Tf. 

On the other hand, with the development of the computer science and the discreteness of 
the micro physics, to look for the integrable discrete forms of the useful continuous integrable 
models becomes a hot topic in nonlinear science. Actually, one continuous integrable model 
may have some different integrable discrete forms. For instance, the integrable discrete 
versions of the KdV equation may be the Volterra equation, the Toda lattice jSj, the Ablowitz 
model 9] and the hybrid lattice etc. 



In section II, we propose a coupled Volterra system which is a discrete version of the 
coupled KdV system ([6])-([7]) and study the integrabihty of the coupled Volterra system. 
The periodic cnoidal waves, solitons (negatons), positons and complexitons are investigated 
by using some suitable solution ansatzs in section III and section IV for a > and a < 0, 
respectively. A short summary is presented in the last section. 

II. A COUPLED VOLTERRA SYSTEM 

It is interesting that the following coupled Volterra system 

ant - «n(an-i - On+l) " - = 0, (8) 

Kt - an{bn~l - K+l) - bn{an~l - «n+l) = 0, (9) 

is an integrable extension of the usual Volterra system ([T]). It is obvious that both 6„ = and 
bn = -^a{n, t) reduce the coupled Volterra system ([8])-([9]) to the usual Volterra equation 

©• 

The integrabihty of the coupled Volterra system ([8])-(l9]) is guaranteed by the following 
theorem. 

Theorem. The coupled Volterra system dH])-© possesses the following Lax pair, 

Lijn = Al/Jn, (10) 
ijnt = Mij^, (11) 



bnT+ anT+ + T_ 
a„a„+i + abnhn+i a(a„6„+i + fenOn+i) 






' Ai a\2 1 








I A2 Ai j 




[ ^2n / 



(12) 

Tl, (13) 



and and T_ are shift operators defined by 

T+fn = fn+U T.fn = fn-1 (14) 

for an arbitrary function fn- 

Proof. By the direct calculations, one can prove that the matrices M and A are commutable 



and then the compatibihty condition of ffT2l) and (fT3|) is just the zero curvature condition 

Lt + LM - ML = 0. (15) 

Substituting the definition equations of L and M into ( fTSll just leads to the coupled Volterra 
system (IHI)-®. The theorem is proved. 

Another interesting fact is that the coupled Volterra system (IHl)^® is really a discrete 
form of the coupled KdV system (I6])-©. Applying the following continuous limiting proce- 
dure 

an = 1 + 6\ (^{n - 2t)S, + 0{S^) 

= l + 6\{x, t) + 0{6''), (16) 
bn = S^v (^{n - 2t)5, ^5h^ + 0{S^) 

= S^ix, T) + OiS^), (17) 
a„±i = l + 6\{x±S, t)+0{6^), (18) 



b^^^ = 6^v{x±5, t)+0{S^) (19) 



to (I8])-(I9]), we have 



^{ur + 6avv^ + 6uu^ + u^^^)S^ + 0{6^) = 0, (20) 
^{vr + Qvu^ + 6uv^ + v^^^)6^ + 0{6^) = (21) 
which is just the special coupled KdV system (I6])-(I7]) with t t. 

III. CNOIDAL WAVES, SOLITONS AND POSITONS OF THE COUPLED 
VOLTERRA SYSTEM FOR a > 

In the continuous case, the rational expansion of the elliptic and hyperbolic functions is 
one of the simplest methods to find travelling periodic and solitary wave solutions. Fortu- 
nately, this method is valid also in discrete cases. 

To find some types of cnoidal wave solutions of the coupled Volterra system ([HJ-Q, one 
may take some types of the elliptic function expansion ansatzs, say, 

ain t) = ^^^'^'^^ (22) 
"^"''^ Rif,g,hy ^''^ 

,( ,x Q{f,g,h) . . 



where {P{f,g,h), Q{f,g,h),R{f,g,h)} are polynomial functions of {/, g, h} and / = 
sn{kn + ct,m), g = cn{kn + ct, m) and h = dn(A;n + ct, m) are Jacobi elliptic functions with 
three constant parameters, the wave number k, the angular frequency u and the modulus 
m. Because of the computational difficulty, here, we just give some special examples of 
(I22])-(I231). 

Case 1. The first simple expansion ansatz reads 

. , fli + Cicn^ffcn + ct, m) . . 

a(n,t) = — -, (24) 

^2 + C2Cn^(K?7, + ct, mj 

A[l + aocn^ikn + ct, m)] 

bin,t) = — — — ^ (25) 

a2 + C2cn"'(A:n + ct, mj 

where ai, ci, a2, C2, A, k, m and c are constants that should be determined later. 

Substituting fl24|) - fl25l) into ([8])-([9]) and vanishing the coefficients of the different powers of 
the Jacobi elliptic function cn^(A;?7, + ct), we can obtain a complicated determining equation 
system for the undetermined constants. Fortunately, there exist a unique general solution 
for the determining equation system. The result reads (s = sn(A;5, m), d = dn{k6, m), C = 
cn{kS, m)) 

2 _ 4a2C2(a2 + C2)[m^(a2 + Ca) - C2] , . 

' ~ K(a2 + C2)2-C2]2 ' ^ > 

.2 _ Q2C2C^(a2 + C2)c^[m^(a2 + C2) - C2] , , 

a(a2«o — C2Y 



ai 



cm^s{a2 + C2){a2 - C2 + 2aoa2)[(a2 + C2)m^s^ - 2C2] ca2 



AdCc2{aoa2 — C2) 2dCs 
^ cc2s{2aoa2 - C2)(2a2 + 2c2 - C2g^) ^^g-j 
4(tC(a2 + C2)(aoa2 - C2) 
^ _ Qi(qoC2 + 2c2 - aoa2) _^ c(a2 + C2 - C2.s^)^(c2 - apaa) 

a2 + 2aoa2 - C2 2dCs{a2 + C2)(a2 + 2aoa2 - C2) ' 

while the constants 02, C2, c, oq and m remain to be free parameters. 

From fl26l) and fl271) . it is known that the real condition of the periodic solution requires 

a > 0. 

Fig. 1 shows the structure of the periodic solution flM|) with the parameter selections 

a = c = a2 = C2 = 1, ao = 2, 5 = l,m = 0.999 (30) 



at time t = 0. 
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FIG. 1: The structure of the periodic wave expressed by ()24p with the parameters (I30p at time 
t = 0. 

It is remarkable that for the a > case, in addition to the above cnoidal wave solution, we 
can find many other nonequivalent periodic waves by similar procedures. In the following, 
we just list some of them. 
Case 2. 

ai + Cicn(kn + ct, m) 

a{n,t) = ^- r, (31) 

0-2 + C2cn[Kn + ct, m) 

Ai + A2CYi{kn + ct, m)] 

o[n,t) = — (321 

a2 + C2cn[kn + ct, m) 

where the constants Oi, Ci, a2, C2, Ai, A2, k, c and m satisfy the following conditions: 

tti = -^^{rn^at + 1 — m?), (33) 



A2 2a2\c + 2cisd — cm? [l 



al 



(34) 



Case 3. 



Al cs^ [m2(l - a2)2 - l] + 2al[c + 2cisrf) ' 

cs2rm2(l-a2)2-ll +2a2(c + 2cisd) , , 

4saa2V« 

C = 1 + [m2(l - a2)2 - 1] , C2 = 1. (36) 
/a2 



a{n,t) = -4^ ■ (37) 

a2 + C2dn(A;?7, + ct, m) 

6(n t) = ^1 + ^2cin(fcn + ct, m)] 
a2 + C2dn{kn + ct, m) 



with the constant constraints 



Case 4. 



SC 

ai = - — —{m^ + a2 - 1) + Cia2, 

A2 2a2[c + 2cisC - cs^{l - al)] 

A " cs2 [(1 - - m2] + 2ai(c + 2cisC) ' 
^ _^cs2[(l-ai)2-m2] +2ai(c + 2cisC) 
^ 4sCa2-\/a 

^ = 1 + ^ [(1 - «2)' - , C2 = 1. 

za2 



, , ai + cisn(A;n + ct, m) 
a{n,t) — 



b{n,t) 



02 + C2sn(kn + ct, m) ' 
Ai + A2sn(kn + ct, m)] 



a2 + C2sn(kn + ct, m) 
with the constant constraints 

ai c(Cd-l) 2 4 IX 

02 2s(l + m^a|) 

A2 _ ( 1 + m^al) (2sai + cdCa2) 

Ai 02 [m?a2{ca2 + 2sai) + 2sai + 2cdCa2 — 002] ' 

m'^a\{ca2 + 2sa\) + 2sai + 2cdCa2 — ca2 
Ai — ± — 



2a2y/2a{l - dC){l 
2_2 al{l-dC ) 
1 + rri^ai 



2 ■""'2V-^ J 1 

s = — ^"1-, C2 = 1. 



Case 5. 



, , aicnf/cn + ct, m) + Cinnikn + ct, m 
a[n, t) — 



b{n, t) 



a2cn{kn + ct, m) + C2sn(A;n + ct, m) ' 
Aicn(/cn + ct, m) + A2sa.{kn + ct, m) 
a2ca{kn + ct, m) + C2sn{kn + ct, m) 

with the constant constraints 

ai c((i — 1) (771^02 — ^2 + 1) 
" ^ ~ 2sC[m2a| - (1 + a2)2] ' 

A2 _ [(1 + a^)2 - m'^aj]{cda2 + 2sCai) 

Ai a2\m?a2{ca2 + 2saiC) — (1 + a2){2sCai + 2c(ia2 — ca2 + ca2 + 2sCaio|)] ' 
, 09(002 + 2saiC) — (1 + a^){2sCai + 2c(ia2 — ca2 + ca, + 2sCaiao) 



Case 6. 



a(n, t) - 
b{n,t) -- 

with the constant constraints 



aidn(A;T2, + ct, m) + Cisn{kn + ct, m) 
a2dn(fcT2, + ct, m) + C2sn(/cn + ct, m) ' 
Aidn{kn + ct, m) + A2Sii{kn + ct, m) 
a2dn(A;n + ct, m) + C2sn(A;?7, + ct, m) 



(55) 
(56) 



Cl 



^9 



ai 

^2 



c(C — \){m'^a\ — 1 — m^flg) 
2s(i[(l + m^Cg)^ — m^ttg] 



[(1 + rrra\) — ■rrral]{cCa2 + 2sdai 



Ai a2 [2sdai[{l + w?alY - m'^al] + ca2[m'^al + m2ai(2C - a^) + 2C - 1 
s[2sdai[{l + m^a^)^ - m^af] + coaim^a^ + m2a2(2C - af) + 2C - 1]] 

2al{l-C) 



(1 + a2m- 



2^2 



9 4 ' 



C2 = l. 



(57) 

(58) 

(59) 
(60) 



Case 7. 



a(n, t) : 
b{n,t) -- 

with the constant constraints 



aidia{kn + ct, m) + ciCYi{kn + ct, m) 
a2dn(A;n + ct, m) + C2ca{kn + ct, m) ' 
Ai(hi{kn + ct, m) + yl2cn(A;n + ct, m) 
a2dn(/i;?7, + ct, m) + C2cn{kn + ct, m) 



c(Cd- l)(m2a| - 1) 
2s(l + m^a^) ' 
A2 _ (1 + m^a^) (cC(ta2 + 2sai) 

v4i 02 \rin?a2{ca2 + 2sai) + 2sai + 2cda2C — ca2] ' 
m^a2(ca2 + 2sai) + 2sai + 2cda2C — ca2 



(61) 
(62) 



Cl = — 
02 



2a2v/2a(l -dC)(l 



2al{l-dC) 
1 + anm^ 



C2 = 1. 



(63) 
(64) 
(65) 
(66) 



Especially, when m ■ 
solutions. For instance, 



1, the previous periodic wave solutions become negaton (soliton) 



and fl25l) become the soliton solution 
Oi + Cisech^(A;n + ct) 



a(n, t) 
6(n,t) 



0'2 + C2sech^(A;n + ct) ' 
74[1 + aosech^(fc?7, + ct)] 
02 + C2sech^(/c?7, + ct) 



(67) 
(68) 
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FIG. 2: The structure of the solitary wave expressed by ()67|) which is the limit case of the figure 
1 for 171 = 1. 



with 



^2 ^ 4c2(a2 + C2) ^2 ^ C2alc^{a2 + C2) 
(a2 + 2c2)2 ' a{a2ao - 02)^ 

_ 2ca2C2(a2 + C2) ajc , . 

S'(a2ao - C2)(a2 + 2C2) 2S'(a2ao - £2)^ ' 



_ 2cc2(a2 + C2) C2c(5a2 + 12a2C2 + 8cl) 

S{a2ao - C2){a2 + 2C2) 2S{a2aQ-C2Y 



Fig. 2 shows the structure of the sohton solution with the same parameter selections 



as (l30l) except for m = 1. 

In [Tj, for the a < case, it is found that the coupled KdV system ([6]) and ([7j) possesses 
not only the analytic solitons (negatons) but also analytical positons and complexitons. Now 
the natural question is whether the coupled Volterra system possesses analytical positons 
and/or complexitons. 

Similarly, all the cnoidal wave solutions presented above will reduced to the positon 
solutions when we take m = 0. For instance, the first type of the positon solutions can be 
obtained from the periodic solution fl2^ and (!25|) by taking m = 0, which has the form of 

_ 2[ai + ci cos^(fcn + ct)] 

' cos{2kn + 2ct) ± cos(M) ' ^ ^ 

^ 2A[l + ao cos\kn + ct)] 

^ v^[cos(2A;ri + 2ct) ± cos(M)] ' ^ 




FIG. 3: A typical positon structure expressed by (f72]l "(f76 ]) with the parameter selections (f77|) . 



where 



ai 



Cl 



[4 + ao(2 + cos{2k5) + cos(4A;5)) T (1 + cto) cos(M) cos(2A;5)] 

4 sin(2/i;5) [2 + oq =F Oo cos(A;(5)] 
c[2 + ao =F flo cos(A;5) cos(2fc5)] 



sin(2/i;5) [2 + oq =F Oo cos(A;5)] 



(74) 
(75) 
(76) 



[2 + ao =F flo cos(A;(5)]^ 
and c, ao and k are arbitrary constants. 

Obviously, the positon solution (!72l) -(173l) is always singular. Fig. 3 shows a special 
structure of (!72|) and (173|) with the parameter selections 



,5 = 1, ao = 2, k = 0.6435011088. 



(77) 



Actually, all the positon solutions obtained from Cases 1 to 7 by taking m = are singular 
except for the trivial constant solution in Case 3. 

For a coupled nonlinear system, there may be different types of solitons and positons. 
The first type of soliton solutions given by ( !67|) and ( !68l) possesses the property that the 
fields a(?7., t) and h{n, t) both have the ring or bell shape. The coupled Volterra system can 
have other kinds of soliton solutions. For instance, by substituting the following solution 
ansatz 

ao + ai cos\i{kn + ct) + a2 cosh(2A;?7, + 2ct) 



a{n, t) 
b{n,t) 



bo + bi cosh{kn + ct) + 62 cosh(2/c?7- + 2ct) ' 

do sinh(/cn + ct) 
bo + bi cosh(/cn + ct) + 62 cosh(2/cn + 2ct) 



(78) 
(79) 



Fig.4a 



Fig.4b 




FIG. 4: Second type of soliton structure expressed by ([78l) - (f85]l with the parameter selections (f86]l . 



into the coupled Volterra system ([8]) and ([9]), one can find that 

do 
bo 



cos -k6 
2 



=Fcsin(co 
^sin(A;5)(l + cos(2co) + cos{k6)), 



cos -kS 
2 



bi = ±4v^sin(A;5) cos(co) cos ^-/c^^ 



ao = -2 



cos(2co) + 2 cos ( 2^*^] '^^^ \2^^ 



ai = =Fccos(co) 
62 = \/asin(A;5), 02 



cos ( -k6] + cos ( -k6 
2 \2 



c 

2' 



(80) 
(81) 
(82) 

(83) 

(84) 
(85) 



where c, Cq and k are arbitrary constants. 

From (!80l) -(l85l). we know that this kind of sohton solution is also valid only for a > 0. 
The solution (!78l) - (179|) may be singular or analytical based on the different selections of the 
parameters. 

Fig. 4 exhibits the structure of the soliton solution given by (!78l) - (!85|) with the special 
parameter selections 



6 = Co = a = 1, c = k 



(86) 



for the upper sign. It is clear that in this case the soliton structures are different for the 
fields a{n, t) and b{n, t). The soliton structure for a{n, t) is still a bell or ring shape while 
that for b{n, t) becomes staggered. 

It is noted that because of the arbitrariness of the constants k, c and Cq, if we take 



then the sohton solution (178!) - (!85l) is transformed to another type of positon solutions 

ttQ + ai cos{kn + ct) + a2 cos{2kn + 2ct) 



a{n, t) 
b{n,t) 



bo + bi cos{kn + ct) + 62 cos(2A;?7, + 2ct) 

do sm{kn + ct) 
bo + bi cos{kn + ct) + 62 cos(2A;n + 2ct) 



(88) 
(89) 



with 



do = =Fcsinh(co) 



cosh ( -k6 I — cosh ( -k6 



bo = A/a sinh(A;5)[l + cosh(2co) + cosh(A;5)], 
bi = ±4^/a smh.{k6) cosh.{co) cosh (-k6 



«o = -< 



cosh(2co) + 2 cosh ( -k6 ) cosh 



-k6 



ai = =Fccosh(co) 



cosh I I ~^ cosh 



-k6 



f'2 = y/a smh{k6) , 02 = —-- 



(90) 
(91) 
(92) 

(93) 

(94) 
(95) 



Fig. 5 displays the structure of the positon solution fl88l) - fl95l) for the upper sign with the 
special parameter selections 



a 



1, c = k 



IT 



-,co 



3' 



(96) 



Evidently, in this case the positon structure is still singular. Up to now, for the a > case, 
we have not yet found any analytical positons and complexitons. However, similar to the 
single KdV equation and the Toda system , we believe that the coupled Volterra system 
dH]) and Q for a > does not possess analytical positons and complexitons. 



IV. SOLITONS, POSITONS AND COMPLEXITONS OF THE COUPLED 
VOLTERRA SYSTEM FOR a < 

For the coupled real Volterra system ([8]) and Qj with a < 0, all the cnoidal wave solution 
forms in the last section are not valid, however, the function expansion ansatz fl78p and (17^ 
is still applicable to obtain some soliton solutions. 

Similar to the last section, after substituting (JTHl) and (1791) into (IHll-Q for a < and 




FIG. 5: A typical positon structure expressed by (f88]l -(f95 ]) with the parameter selections ([96 



solving the determining equations of the parameters, one can find that 

2csin(co) sinh {^k6) 



-a 



dn 



61 = 4cos(co) cosh ( -k5 ) , 



1 + cos(2co) + cosh(A;5), 



smh.{k6) 
ccos(co) 



[cos(2co) + cosh{kS) + cosh(2M)], 



smh.{kS) 
62 = 1, 02 



c 



(97) 
(98) 
(99) 

(100) 

(101) 
(102) 



2 sinh(A;5) ' 

where c, Cq and k are arbitrary constants. 

Fig. 6 shows the structure of the soliton solution expressed by (|78l) and (1791) with (1971) - 
( 11021) and the parameter selections 



TV , 1 e 

Co = — , c = k = —, = 1, a = —1. 



(103) 



In the continuous case, it has been proved that the coupled KdV system ([6]) and ([7]) with 
a < have some types of analytical positon and complexiton solutions 7|]. It is interesting 
that in the discrete case, the coupled Volterra system also possesses analytical positons and 
complexitons. 

To get analytical positon solutions, one can directly apply the constant transformations 



a(n,t) 



-1.02 



Fig.6a 



Fig.6b 




FIG. 6: The structures of the soliton expressed by (a) ([78]) and (b) ([79j) with (f971)- (fT02]) and the 
parameter selections ()103p . 



to ([78D-([79D with (|97|)- ffT02|) . The result still has the form ([88])-([89]) but with the constants 

2csinh(co) sin 



-a 



bo = 1 + cosh(2co) + cos(fc5), 



bi = 4 cosh(co) cos ( -kd 



ai 
b2 



2 sm{k6) 
ccosh(co) 
sm{k6) 

1, 02 -- 



cosh(2co) + cos{k6) + cos{2k6)], 



cos ( 2^*^] '"^^ \2^^ 



(104) 
(105) 
(106) 

(107) 

(108) 
(109) 



2sin(M)' 

Fig. 7 reveals the structure of the analytical positon solution expressed by (ITHjl and (179|) 
with fll04l) - fll09p and the parameter selections 

Co = 1, c = 1, k = 2, 6=1, a = -l. (110) 

To get analytical complexiton solutions of the coupled Volterra system ([S])-(IH]), we may 
use the following solution ansatz 

_ gp + ai cosh(gi) cos(^2) + ^2 sinh(gi) sin(^2) + a3[cosh(2gi) + cos(2g2)] 
^ ~ bo + b, cosh(ei) cos(6) + &2 sinh(ei) sin(e2) + 63[cosh(2ei) + cos(26)] ' ^ ^ 
_ rfp + di cosh(^i) cos(g2) + 4 sinh(^i) sin(^2) + (i3[cosh(2^2) + cos(2gi)] 



bo + bi cosh(^i) cos(6) + &2 sinh(^i) sin(^2) + b3[cosh(2^i) + cos(2^2)] 



where 



= kiU + Cit + ^oi, ^ = 1, 2. 



(113) 
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FIG. 7: The structure of the positon expressed by (a) ([78]) and (b) ([79]) with (fTnD - (fT09D and the 
parameter selections (jllOp . 

After finishing tedious calculations, we find that ([8])- ([9]) really possesses analytical complexi- 
ton solutions (an analytical complexiton is just a usual breather) if a < while the constants 
ttj, bj and Cj for j = 0, 1, 2 and 3 should be determined by 

flQ = Ci[cos(4c) sinh(26) — sinh(66) — cos(2c) sinh(46)] 

+C2[sin(2c) cosh(46) - sin(4c) cosh(26) - sin(6c)], (114) 
ai = =f2 {[sinh(56) cos(c) — cos(5c) sinh(6) + cos(c) sinh(6) + sinh(36) cos(3c)]ci 

+ [sin(c) cosh(b) — sin(c) cosh(56) + sin(5c) cosh(6) + sin(3c) cosh(36)]c2} , (115) 
02 = =f2 {[sin(3c) cosh(36) — sin(5c) cosh(6) + sin(c) cosh(6) + sin(c) cosh(56)]ci 

-|-[sinh(56) cos(c) — sinh(36) cos(3c) — cos(5c) sinh(6) — cos(c) sinh(6)]c2} , (116) 



as = 


— 2c2 cosh(26) sin(2c) — 2ci sinh(26) cos(2c), 




(117) 


bo = - 


1 

^ [cos(2c)(l + 2cosh(46)) + cosh(26)(l 2cos(4c)) + 

\/—a 


cosh (66) 


-cos(6c)](118) 


bi = - 


=F4 

—j= {cosh(6) cos(5c) + cosh(6) cos(3c) — cos(c) cosh(56) 
\/—a 


— cos(c) 


cosh(3fe)},(119) 


b2 = - 


±4 

—^= {— sinh(6) sin(5c) + sinh(6) sin(3c) + sin(c) sinh(56) 
\/—cy 


— sin(c) 


sinh(36)},(120) 


bs = - 


2 

-^=[cosh(46) - cos(4c)], 




(121) 




FIG. 8: The structure of the complexiton expressed by (a) (fTTTI) and (b) (fTT2D with (fTTil) - (fT25]l 
and the parameter selections ()126p . 

do = — C2 cos(4c) sinh(25) — ci sin(4c) cosh(2fe) — ci sin(6c) 

+C2 sinh(66) + Ci sin(2c) cosh(46) + Ca cos(2c) sinh(4fe), (122) 
di = ±2 {[sinh(56) cos(c) — cos(5c) sinh(6) + cos(c) sinh(6) + sinh(36) cos(3c)]ci 

+ [sin(c) cosh(6) — sin(c) cosh(56) + sin(5c) cosh(6) + sin(3c) cosh(36)]c2} , (123) 
d2 = ±2 {[cos(5c) sinh(6) + cos(c) sinh(6) — sinh(56) cos(c) + sinh(36) cos(3c)]ci 

-|-[sin(3c) cosh(36) — sin(5c) cosh(6) + sin(c) cosh(6) + sin(c) cosh(56)]c2} , (124) 
ds = 2c2 cos(2c) sinh(26) - 2ci sin(2c) cosh(26), (125) 

with 

b = ^ki6, c = ^k26 

and fci, ^21 Ci, ^01 and ^02 being arbitrary constants. 

Fig. 8 displays the complexiton structure expressed by fillip and (11121) with flll4p -( fT25ll 
and the special parameter selections 

k2 = ^01 = 0, h = 0.1, ci = 0.2, C2 = 0.7, eo2 = ^,5 = 1. (126) 

V. SUMMARY AND DISCUSSIONS 

In summary, the special coupled integrable KdV system (l6])-([7j) is discreterized to an 
integrable coupled Volterra system. The Lax integrability of the coupled Volterra system 
is proved. By using a simple rational expansion method of the Jacobi elliptic functions, 
trigonometric functions and hyperbolic functions, various exact solutions are found. 



For the coupled Vorterra system ([HJ-Q with a > there are many types of cnoidal 
waves described by different types of Jacobi elhptic functions. On one hand, whence the 
modulus of the cnoidal wave tends to 1, the wave tends to a negaton solution. If a negaton 
is analytical, then it is called soliton (or solitary wave for nonintegrable systems). It is found 
that only one of these cnoidal waves can be reduced to a single analytical negaton solution 
when the modulus, m, of the model tends to 1. The soliton solution has a ring or bell shape 
for both fields a{n,t) and b{n,t). On the other hand, whence the modulus of the cnoidal 
wave tends to 0, the wave tends to a positon solution. It is found that all the nontrivial 
positons obtained from the cnoidal waves of section III are singular. This fact is similar to 

rn 

the KdV and Toda system 

In the a > case, there are two types of single soliton solutions. In addition to the 
above mentioned analytical negaton solution, there is a different type of solitons which has 
different shapes for the fields a{n, t) and b{n, t). The new type of soliton solution is obtained 
by taking a more complicated rational expansion of the hyperbolic functions. Because of 
the calculation difficulty, we have not yet found any cnoidal wave extension of this type 
of negaton solution even utilizing computer algebras. The field a{n, t) for the second type 
of soliton solution also possesses the ring or bell shape while the field b{n, t) possesses a 
staggered shape. 

For the coupled Vorterra system ([HJ-Q with a < 0, though we have not yet found the 
cnoidal wave solutions, many other kinds of physically significant solutions, such as the 
solitons, analytical positons and analytical complexitons are found. The structure of the 
soliton solution in this case is similar to that of the second type of the solitons for a > 0. 

It is also interesting that the positon solution can be obtained by many methods. In this 
paper, we demonstrate that the analytical positon solutions of the coupled Vorterra system 
©-([I]) with a < can be simply obtained from the negatons by means of the constant 
analytical extensions. 

To obtain analytical complexiton solutions of the coupled Vorterra system (JHl)^© with 
a < 0, a more complicated rational expansion of both the hyperbolic functions and trigono- 
metric functions is used. 

Finally it is worth to indicated that the real solutions of the coupled Vorterra system ([8])- 
([9]) with a < can also be obtained by means of the analytical continuous extensions and vice 
versa. Some analytical continuous extension examples have been given in section IV. Here, 
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FIG. 9: The evolution of the two-sohton interaction solution expressed by (a) and (b) ()112|] 

with the transformations (jl27p and the parameter selections (|128p . 

we just mention a further interesting example. If we apply the constant transformations 

h V^h, C2 V^C2, ^02 V^Co2, (127) 



to flllip - (lll2p . then the complexiton solution of a < case becomes a two-soliton solution 
for a > 0. 

Fig. 9 shows the two-soliton interaction related to flllip - flll2p with the transformation 
(11270 under the parameter (after (11270 ) selections 



TT 



h = 6i = 0, ki = 0.1, ci = 0.1, C2 = 1,^02 = ^,5 = 1. 



(128) 



Though we have obtained many types of exact solutions of the model via a simple function 
expansion method, various problems, such as the general multiple soliton solutions and r 
function solutions are still open. As a discrete form of the significant physical model, the 
more about the model will be studied further. 
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